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ABSTRACT 

We report the discovery of two concentric partial Einstein rings around the gravitational lens 
SDSSJ0946+1006, as part of the Sloan Lens ACS Survey. The main lens is at redshift z\ = 0.222, 
while the inner ring (1) is at redshift z s \ — 0.609 and Einstein radius i?Eini = 1-43 ± 0.01". The 
wider image separation (i?Ein2 — 2.07±0.02") of the outer ring (2) implies that it is at higher redshift 
than Ring 1. Although no spectroscopic feature was detected in ~ 9 hours of spectroscopy at the 
Keck I Telescope, the detection of Ring 2 in the F814W ACS filter implies an upper limit on the 
redshift of z S 2 < 6.9. The lens configuration can be well described by a power law total mass density 
profile for the main lens ptot oc r 1 with logarithmic slope 7' = 2.00 ± 0.03 (i.e. close to isothermal), 
velocity dispersion ctsie = 287 ± 5kms~ (in good agreement with the stellar velocity dispersion 
cr,,,* = 284 ± 24kms~ 1 ) with little dependence upon cosmological parameters or the redshift of Ring 
2. Using strong lensing constraints only we show that the enclosed mass to light ratio increases as a 
function of radius, inconsistent with mass following light. Adopting a prior on the stellar mass to light 
ratio from previous SLACS work we infer that 73 ± 9% of the mass is in form of dark matter within 
the cylinder of radius equal to the effective radius of the lens. We consider whether the double source 
plane configuration can be used to constrain cosmological parameters exploiting the ratios of angular 
distance ratios entering the set of lens equations. We find that constraints for SDSSJ0946+1006 are 
uninteresting due to the sub-optimal lens and source redshifts for this application. We then consider 
the perturbing effect of the mass associated with Ring 1 (modeled as a singular isothermal sphere) 
building a double lens plane compound lens model. This introduces minor changes to the mass of the 
main lens, allows to estimate the redshift of the Ring 2 (z S 2 = 3.1^'q), and the mass of the source 

responsible for Ring 1 (osie.sI = 941^ kms^ 1 ). We conclude by examining the prospects of doing 
cosmography with a sample of 50 double source plane gravitational lenses, expected from future space 
based surveys such as DUNE or JDEM. Taking full account of the uncertainty in the mass density 
profile of the main lens, and of the effect of the perturber, and assuming known redshifts for both 
sources, we find that such a sample could be used to measure f2 m and w with 10% accuracy, assuming 
a flat cosmological model. 

Subject headings: Gravitational lensing - galaxies : Ellipticals and lenticulars, cD - galaxies: structure 
- galaxies: halos - cosmology: dark matter - cosmology: cosmological parameters 
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Measuring the mass distribution of galaxies is essential 
for understanding a variety of astrophysical processes. 
Extended mass profiles of galaxies provide evidence for 
dark matter either using rotation curves (e.g. Rubin et al. 
1980; van Albada et al. 1985; Swaters et al. 2003), weak 
lensing (e.g. Brainerd et al. 1996; Hoekstra et al. 2004; 
Sheldon et al. 2004; Mandelbaum et al. 2006), or dy- 
namics of satellite galaxies (e.g. Prada et al. 2003; Con- 
roy et al. 2007) which is one of the main ingredients of 
the standard A cold dark matter (ACDM) cosmologi- 
cal model. At galactic and subgalactic scales, numerical 
cosmological simulations make quantitative predictions 
regarding, e.g., the inner slope of mass density profiles 
and the existence of dark matter substructure. Precise 
mass measurements are key to test the predictions and 
provide empirical input to further improve the models. 

Gravitational lensing has emerged in the last two 
decades as one of the most powerful ways to measure 
the mass distributions of galaxies, by itself or in com- 
bination with other diagnostics. Although strong grav- 
itational lenses are relatively rare in the sky (< 20 per 
square degree at space-based depth and resolution; Mar- 
shall et al. 2005; Moustakas et al. 2007), the number 
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of known galaxy-scale gravitational lens systems has in- 
creased well beyond a hundred as a result of a number of 
dedicated efforts exploiting a variety of techniques (e.g. 
Warren et al. 1996; Ratnatunga et al. 1999; Kochanek 
et al. 1999; Myers et al. 2003; Bolton et al. 2004; Ca- 
banac et al. 2007). The increased number of systems, 
together with the improvement of modeling techniques 
(e.g. Kochanek & Narayan 1992; Warren & Dye 2003; 
Treu & Koopmans 2004; Brewer & Lewis 2006; Suyu 
et al. 2006; Wayth & Webster 2006; Barnabe & Koop- 
mans 2007), has not only enabled considerable progress 
in the use of this diagnostic for the study of the mass 
distribution of early and most recently late-type galax- 
ies, but also for cosmography, i.e. the determination of 
cosmological parameters (e.g. Golse et al. 2002; Soucail 
et al. 2004; Dalai et al. 2005). 

Given the already small optical depth for strong lens- 
ing, the lensing of multiple background sources by a 
single foreground galaxy is an extremely rare event. 
At Hubble Space Telescope (HST) resolution (FWHM 
<~ 0'.'12) and depth (Lab ~ 27) it is expected that one 
massive early-type galaxy (which dominate the lensing 
cross-section) in about 200 is a strong lens (Marshall 
et al. 2005). Taking into account the strong dependence 
of the lensing cross-section on lens galaxy velocity dis- 
persion (cx (T 4 ), and the population of lens galaxies, we 
estimate that about one lens galaxy in ~ 40 — 80 could 
be a double source plane strong gravitational lens (see 
appendix A). For these reasons, at most a handful of 
double lenses are to be found in the largest spectroscopic 
surveys of early-type galaxies such as the luminous red 
galaxies of the Sloan Digital Sky Survey. However, future 
high resolution imaging surveys such as those planned for 
JDEM and DUNE (Aldering & the SNAP collaboration 
2004; Refregier et al. 2006) will increase the number of 
known lenses by 2-3 orders of magnitude (Marshall et al. 
2005), and hence should be able to provide large statis- 
tical samples of double source plane gravitational lenses, 
opening up the possibility of qualitatively new applica- 
tions of gravitational lensing for the study of galaxy for- 
mation and cosmography. 

We report here the discovery of the first double source 
plane partial Einstein Ring. The gravitational lens sys- 
tem SDSSJ0946+1006 , was discovered as part of the 
Sloan Lens ACS (SLACS) Survey (Bolton et al. 2005, 
2006a; Treu et al. 2006; Koopmans et al. 2006; Bolton 
et al. 2007; Gavazzi et al. 2007). The object was first se- 
lected by the presence of multiple emission lines at higher 
redshift in the residuals of an absorption line spectrum 
from the SDSS database as described by Bolton et al. 
(2004) and then confirmed as a strong lens by high res- 
olution imaging with the Advanced Camera for Surveys 
aboard HST. In addition to an Einstein ring due to the 
source (hereafter source 1) responsible for the emission 
lines detected in the SDSS spectrum, the Hubble image 
also shows a second multiply imaged system forming a 
broken Einstein Ring with a larger diameter then the in- 
ner ring (hereafter source 2). This configuration can only 
arise if the two lensed systems are at different redshifts 
and well aligned with the center of the lensing galaxy. It 
is a great opportunity that a double source plane lens has 
been found among the approximately 90 lenses discov- 
ered by the SLACS collaboration to date (Bolton et al. 
2008, in prep.). 



The goal of this paper is to study and model this pe- 
culiar system in detail, as an illustration of some astro- 
physical applications of double source plane compound 
lenses, including i) the determination of the mass density 
profile of the lens galaxy independent of dynamical con- 
straints; ii) placing limits on the mass of source 1 based 
on multiple lens plane modeling; iii) estimating the red- 
shift of source 2 and the cosmological parameters from 
the angular distance size ratios. The paper is therefore 
organized as follows. Section 2 summarizes the observa- 
tions, photometric and spectroscopic measurements, and 
discusses the morphology of the lens system. Section 
3 describes our gravitational lens modeling methodol- 
ogy. Section 4 gives the main results in terms of con- 
straints on the mass distribution of the lens galaxy and 
of source 1. Section 5 discusses the use of double source 
plane lenses as a tool for cosmography using the example 
of SDSSJ0946+1006 and also addresses the potential of 
large samples of such double source plane lenses for the 
same purpose. In section 6 we summarize our results and 
briefly conclude. 

Unless otherwise stated we assume a concordance cos- 
mology with H = 70/i7o kms _1 Mpc _1 , £! m = 0.3 and 
£1a = 0.7. All magnitudes are expressed in the AB sys- 
tem. 

2. DATA 

The lens galaxy SDSS J0946+ 1006 was first identified 
in the spectroscopic SDSS database based on the redshift 
of the lensing galaxy z\ = 0.222 and that of a background 
source at z s \ = 0.609 (hereafter source 1), as described 
by Bolton et al. (2004, 2006a), and Bolton et al. (2008, 
in prep.). This section describes HST follow-up imaging 
(§2.1), the properties of the lens (§2.2) and lensed (§2.3) 
galaxies. 

2.1. Hubble Space Telescope observations & data 
reduction 

SDSSJ0946+1006 was then imaged with the ACS on 
board the HST (cycle 15, Prog. 10886, PI Bolton). 
The Wide Field Channel with filter F814W was used 
for a total exposure time of 2096 s. Four sub- 
exposures were obtained with a semi-integer pixel offset 
(acs-wf c-dither-box) to ensure proper cosmic ray re- 
moval and sampling of the point spread function. The 
image reduction process is described in (Gavazzi et al. 
2007) and results in a O'.'03/pixel spatial sampling. This 
pixel size provides good sampling of the PSF for weak 
lensing applications, at the (small) price of inducing noise 
correlation over scales of 1-2 pixels. This is accounted for 
in our analysis by correcting pixel variances according to 
the procedure described by Casertano et al. (2000). 

Figure 1 shows the HST image of the lens galaxy field 
together with an enlarged view of the lensed features, 
after subtraction of a smooth model for the lens surface 
brightness distribution. For reference, one arcsecond in 
the lens plane subtends a physical scale of 3.580 /if kpc. 

2.2. Lens galaxy properties 

The two-dimensional lens surface brightness was fit- 
ted with galf it (Peng et al. 2002) using two elliptical 
Sersic components. The addition of a second component 
is needed to provide a good fit in the center, and to repro- 
duce the isophotal twist in the outer regions. To reduce 
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Fig. 1. — HST/F814W overview of the lens system SDSSJ0946+1006. The right panel is a zoom onto the lens showing two 
concentric partial ring-like structures after subtracting the lens surface brightness. 



the effect of lensed features in the fit we proceeded it- 
eratively. We first masked the lensed features manually, 
then we performed galfit fits creating masks by 4-ct 
clipping. Two iterations were needed to achieve conver- 
gence. 

The total magnitude of the lens obtained by sum- 
ming the flux of the two Sersic models is F814W" = 
17.110 ± 0.002 after correction for Galactic extinction 
(Schlcgel et al. 1998). The rest- frame V band abso- 
lute magnitude is M v = -22.286 ± 0.025 using the K- 
correction of Treu et al. (2006). The errors are dom- 
inated by systematic uncertainties on the K-correction 
term. The most concentrated Sersic component C\ dom- 
inates at the center and accounts for about 17.5% of the 
total lens flux. The effective radius of c\ is about 0.4" 
whereas that of C2 is ~ 3" with about 10% relative ac- 
curacy. Similarly, the Sersic indexes are n Cl ~ 1.23 and 
n C2 ~ 1.75. 

To measure the one dimensional light profile of the 
lens galaxy, we used the IRAF task ellipse. Fig. 2 
shows the radial change of ellipticity and position angle 
of the light distribution. There is a clear indication of 
a sharp change in position angle and ellipticity between 
1 — 2". This isophotal twist is well captured by the double 
Sersic profile fit, that requires different PAs for the two 
components. Therefore we conclude that the lens galaxy 
is made of two misaligned components, having similar 
surface brightness at radius ~ 0.6". 

For comparison, a single component Sersic fit yields 
n ~ 3.73, consistent with the typical light profiles of 
massive early-type galaxies. The effective radius of the 
composite surface brightness distribution is found to be 
R cS = 2.02 ± 0.10 arcsec ~ 7.29 ± 0.37 kpc, where 
we assumed a typical relative uncertainty of about 5% 
as discussed in (Treu et al. 2006). It is also consistent 
with an independent measurement reported by Bolton 



et al. (2008, in prep.) who considered de Vaucouleurs 
surface brightness distributions (n = 4 by construction) . 
Note that we use the same convention for all characteris- 
tic radii reported throughout. For elliptical distributions 
radii are expressed at the intermediate radius {i.e. the 
geometric mean radius r = yab). 

In addition, the stellar velocity dispersion cr ap = 263 ± 
21kms _1 was measured with SDSS spectroscopy within 
a 3" diameter fiber. We convert this velocity dispersion 
(Tap into the fiducial velocity cr„,* that enters Fundamen- 
tal Plane analyses and measured in an aperture of size 
i? eff /8 using the relation a v< */a ap = (-R e ff/8/i?ap)" 04 ^ 
1.08 (see Treu et al. 2006, and references therein) 

Based on photometric redshifts available online on the 
SDSS webpage (Oyaizu et al. 2007), we note that the 
lens galaxy is the brightest galaxy in its neighborhood. 
Another bright galaxy about 40 arcsec south-west of 
SDSSJ0946+1006 exhibits perturbed isophotes (an ex- 
tended plume) suggesting that it may have flown by re- 
cently and might end up merging onto the lens galaxy. 
Its photometric redshift is z p hot = 0.20 ± 0.04 consistent 
with SDSSJ0946+1006 redshift. The extended envelope 
captured by the double Sersic component fit also sup- 
ports the recent flyby hypothesis (e.g. Bell et al. 2006). 

2.3. Lensed structures 

Two concentric partial ring-like structures are clearly 
seen at radii 1.43 ±0.01" and 2.07 ±0.02" from the cen- 
ter of the lens galaxy (Figure 1). Such a peculiar lensing 
configuration - with widely different image separations of 
nearly concentric multiple image systems - implies that 
the rings come from two sources at different redshift, the 
innermost (Ring 1) corresponding to the nearest back- 
ground source 1 and the outermost (Ring 2) being sig- 
nificantly further away along the optical axis. 

Ring 1 has a typical cusp configuration with 3 merg- 
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Fig. 2. — Results for isophotal fit with IRAF/ellipse. Top 
panel: Position angle versus radius. Middle panel: Axis 
ratio versus radius. The vertical lines show the location of the 
inner and outer Einstein rings which were masked out during 
the fitting process. We also overlay in the top and middle 
panels as a blue solid line the ellipse output performed on 
the best fit galfit two-dimensional brightness distribution. 
Bottom panel: best fit Sersic profiles obtained with galfit. 
The formal error bars on the surface brightness profile are 
smaller than the data points. 

ing conjugate images and a countcrimage on the opposite 
side of the lens and closer to the center than the large 
cusp "arc". This constrains the orientation of the lens 
potential major axis to pass almost through the middle 
of both arcs. Ring 1 is among the brightest ones to have 
been discovered in the SLACS survey (See Bolton et al. 
2008, in prep., for the latest compilation). The observed 
F814W magnitude is mi = 19. 784±0. 006 (extinction cor- 
rected) . The error bar includes only statistical uncertain- 
ties. An additional systematic error of order < 0.1 mag 
is likely present due to uncertainties in the lens galaxy 
subtraction (Marshall et al. 2007). 

Ring 2 presents a nearly symmetrical Einstein cross 
configuration (with a faint bridge between the north and 
west images), implying that the source must lie very close 
to the optical axis. The observed F814W magnitude is 
ni2 = 23.68 ± 0.09, making it about 36 times fainter 
than Ring 1. As for Ring 1 the error bar includes only 
statistical uncertainties. 

No evidence of Ring 2 is present in the SDSS spectrum. 
This can be explained by the low peak surface bright- 
ness of Ring 2 (~ 23 mag arcsec -2 ) and less importantly 
by the diameter of the second ring being slighty larger 
than the 3" SDSS fiber (although see Bolton et al. 2006b, 
for a successful redshift measurement in a similar case). 
Deeper longslit spectroscopy was obtained at Keck Ob- 
servatory with the Low Resolution Imager Spectrograph 
(LRIS) instrument on December 22-23 2006, the total 
integration time being about 9 hours. The goal was two- 
fold: i) obtain the redshift of Ring 2; ii) measure the 
stellar velocity dispersion profile of the main lens. This 
latter aspect will be presented elsewhere. Despite the 
large integration time, we could not measure the source 
redshift z S 2 due to a lack of emission lines in the range 
[3500, 8600A] that do not belong to Ring 1. Since Ring 



2 is detected in the ACS/F814W filter, we can set an 
upper limit on its redshift z S 2 < 6.9 by requiring that 
the Lyman break be at shorter wavelengths than the red 
cutoff of the filter. 

3. LENS MODELING 

3.1. Model definition 

This section describes our adopted strategy to model 
this exceptional lens system. We begin with a simplifying 
assumption. Although the gravitational potential arises 
from both a stellar and a dark matter component, a single 
power law model for the total density profile turns out to 
be a good description of SLACS lenses (Koopmans et al. 
2006). Therefore, we assume the total convergence for a 
source at redshift z s to be of the form: 



re(r, z s ) 



(x 2 + y 2 /q 2 ) 



(1-7O/2 A s m 

D ' [ ' 



with 4 free parameters: the overall normalization b, the 
logarithmic slope of the density profile 7', the axis ra- 
tio q and position angle PAq (omitted in Eq. (1) for 
simplicity) of iso-K ellipses. The familiar case of the 
singular isothermal sphere is that corresponding to a 
slope 7' = 2 and q = 1. In this case relates 
to the velocity dispersion of the isothermal profile by 
boo = 47t((7sie/c) 2 = (trsiE/186.2kms -1 ) 2 arcsec. Note 
that ctsie is nothing but a way of redefining the normal- 
ization of the convergence profile and does not necessar- 
ily correspond in a straightforward sense to the velocity 
dispersion of stars in the lens galaxy. In general, for ev- 
ery combination of model parameters, the stellar velocity 
dispersion of a specified tracer embedded in the potential 
can be computed by solving the Jeans equation and will 
be a function of radius and observational effects such as 
aperture and seeing. 

No assumptions are made about the orientation of the 
position angle PAq of the lens potential. In addition, 
we allow for external shear with modulus 7 cx t and po- 
sition angle PA ext - For a multiple source plane system, 
it is necessary to define a lens plane from which the ex- 
ternal shear comes from since shear has to be scaled by 
the apropriate D\ s /D os term for each source plane. For 
simplicity we assume that the global effect of external 
pertubations comes from the same lens plane z\ — 0.222. 
We expect a strong degeneracy between internal clliptic- 
ity and external shear but include this extra degree of 
freedom in the model to account for any putative twist 
of isopotentials, as suggested by the observed isophotal 
twist in the lens galaxy surface brightness. Note also 
that the need of being able to handle two distinct source 
planes led us to the somewhat unusual definition of 600 
in Eq. (1). With this convention, {b o\fq) 1 ~ 1 D\ s /D os 
is the quantity closest to the &sie (or i?Einst) parame- 
ter used in other SLACS papers (Koopmans et al. 2006; 
Bolton et al. 2008). Note also that the center of mass is 
assumed to match exactly the lens galaxy center of light. 
The unknown redshift of source 2 is also treated as a free 
parameter, for which we assign a uniform 1 < z S 2 < 6.9 
prior. Altogether, we use 7 free parameters to character- 
ize the potential of SDSSJ0946+1006: b x , 7', q, PA 0l 
7o Xt , PA CKt and z s2 . 

In this section and the next, we neglect the extra fo- 
cusing effect of Ring 1 acting as a lens on Ring 2, leaving 
the discussion of this perturber for Section 5. 
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3.2. Methods 

We consider three strategies for studying gravitational 
lens systems with spatially resolved multiple images. 

The first one consists of identifying conjugate bright 
spots in the multiple images and minimizing the dis- 
tance of conjugate points in the source plane. This ap- 
proach is statistically conservative in the sense that it 
only takes partial advantage of the large amount of in- 
formation present in the deep HST data. However, it has 
the benefit of being robust and relatively insensitive to 
the details of the source morphology, and other concerns 
that affect different alternative techniques in the case. 

The second approach is the linear source inversion and 
parametric potential fitting method described by War- 
ren & Dye (2003), Treu & Koopmans (2004), Koopmans 
(2005) and Suyu et al. (2006). A strong advantage of this 
method is that it takes fully into account the amount 
of information contained in each pixel. Although this 
method is robust, there are many degrees of freedom to 
model the intrinsic source surface brightness distribution 
and thus some form of regularisation is needed to avoid 
fitting the noise as described in the references above. 

The third method (e.g. Marshall et al. 2007; Bolton 
et al. 2007, 2008) describes the source as one or several 
components parameterized with elliptical surface bright- 
ness profiles (usually Sersic). In general, this method 
provides good fits to the data as long as not too many 
such components are needed to represent the source, and 
directly provides physically meaningful parameters for 
the source. For high signal-to-noise ratio images of com- 
plex lensed features the dimensionality of the problem 
may increase very fast. 

In the case of a multiple source plane system, two diffi- 
culties arise when using the second and third techniques. 
1) Our current pixellized method does not handle mul- 
tiple source planes (see e.g. Dye et al. 2007, for recent 
progress along this line). 2) The statistical weight given 
to each of the partial rings depends essentially on their 
relative brightness. Since Ring 1 is 36 times brighter 
than Ring 2, it completely dominates the fit. This has 
the unwanted side effect that a physically uninteresting 
morphological mismatch of the inner ring, due for exam- 
ple to poor modeling of the source or of the point spread 
function, overwhelms any mismatch in the physically im- 
portant image separation of the outer ring. 

The goal of the present analysis is to confirm that 
SDSSJ0946+1006 is the first example of a galaxy-scale 
double source plane system and illustrate what kind of 
information can be inferred from such a configuration. 
After experimenting with all three techniques - and in 
light of the difficulties described above - we decided to 
focus on the more straightforward conjugate points mod- 
eling technique, using the other techniques to aid in our 
modeling. 

In practice, the modeling technique adopted here is 
similar to the one used by Gavazzi et al. (2003). The 
merging cusp nature of ring 1 makes the identification of 
quadruply imaged spots hazardous along the elongated 
arc but identifications are much easier between the oppo- 
site counter-image and the elongated arc. The identifi- 
cation of the brightness peak S2 in Ring 2 is obvious. To 
guide the identification process, we also used fits based 
on the pixellized source inversion. We ended up having 



4 spots identified in Ring 1, two of them having 3 clear 
conjugations (Sla, Sic) whereas the other two have only 
have 2 (Sib, Sid). One single bright spot in Ring 2 is im- 
aged 4 times. The typical rms error made on the location 
of spots estimated to be 0.03". Table 1 summarizes the 
coordinates of matched points in the same frame as Fig. 
1. For each knot Sla, Sib, Sic, Sid and S2, multiple 
images with positive parity have an odd labelling num- 
ber. To guide the fitting procedure we also demand the 
image parity to be preserved by the model. Therefore, 
taking into account the unknown position of these spots 
in the source plane, we end up having 18 constraints (see 
Gavazzi et al. 2003) whereas the considered model has 7 
free parameters. Hence the optimization problem has 11 
degrees of freedom. 



TABLE 1 

Summary of pixel coordinates used for lens modeling. 





Img. 1 


Img. 2 


Img. 3 


Img. 4 


Sla 


0.34, -1.50 


-0.94, 0.68 


1.52, 0.19 




Sib 




-1.16, 0.22 


1.44, 0.88 




Sic 


-0.43, -1.42 


-1.10, 0.67 


1.23, 0.88 




Sid 


-0.14, -1.68 


-0.57, 0.96 






S2 


-1.51, -1.78 


1.56,-1.19 


1.55, 1.65 


-1.34, 1.32 



Positions (x,y) of each multiple knot are expressed in arcsec 
(typical rms error 0.03") relative to the lens galaxy surface 
brightness peak (got from galfit modeling, see §2.2). The 
frame position angle is 161.348° relative to North direction. 



4. MODELING RESULTS 

The optimization process and the exploration of the 
parameter space were performed by sampling the poste- 
rior probability distribution function with Monte-Carlo 
Markov Chains (MCMC). We assumed flat priors. Table 
2 summarizes the results ("best fit" values are defined 
as the median value of the marginalized PDF) and their 
corresponding 68% CL uncertainties after marginalizing 
the posterior over all the other parameters. The best fit 
model yields a x 2 /dof = 13.2/11 ~ 1.20 which is statis- 
tically reasonable 10 . 



TABLE 2 

Best-fit model parameters for SDSSJ0946+1006 using a 
single lens plane. 



[arcsec] 



2.54 ± 0.09 
2.00 ± 0.03 
0.869 
-11.8 
0.067 
-31.5 
5.30 



0.017 
-0.013 
+7.0 
8.9 
0.010 
-0.007 
6.9 
4.8 
+ 1.03 
-1.00 

287.0+53 



axis ratio q 
PA 

7cxt 

PA cxt 

Zs2 

ct sie [kms -1 ] 

"unlensed" apparent F814 s i [mag] 22.76 ± 0.02 ± 0.10 

"unlensed" absolute V s i [mag] -19.79 ± 0.05 ± 0.10 

"unlensed" apparent F814 s2 [mag] 27.01 ± 0.09 ± 0.10 

Best fit model parameters and 68.4% confidence limits. Er- 
rors on magnitudes distinguish statistical (first) and system- 
atic from lens light subtraction (second). Angles are in de- 
grees oriented from North to East. 

10 A x 2 distribution with 11 degrees of freedom gives a proba- 
bility of 28% that the \ 2 value will be greater than 13.2 . 
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Fig. 3. — Best-fit single lens plane model for the lens SDSSJ0946+1006. The model parameters were found using the 
identification of conjugate bright knots but the quality of the model is illustrated with a pixelised source inversion technique. 
Top left: observation with the lens light profile subtracted off. Top middle: model prediction in the image plane and associated 
residuals {Top right). The model also predicts the light distribution in the source planes z Bl and z B 2 (Bottom left and right 
respectively). Note a different color stretching for source plane 2 (factor 6) in this latter case. Critical and caustic lines 
corresponding to the two source planes are overlaid (smaller blue for z B i = 0.609 and wider red for z B 2 = 5). 

The results of the best fit model inferred from the con- 
jugation of bright spots is shown in Fig. 3 where we used 
the pixellized source inversion technique to illustrate the 
quality of the fit and the reliability of the conjugation 
method. Although the surface brightness of Ring 1 and 
Ring 2 identified by separate annuli in the image plane 
are inverted separately model predictions in the image 
plane are recombined for convenience. The two source 
planes z sl = 0.609 and z S 2 ~ 5 are also shown. 

As expected, there is a degeneracy between external 
shear and ellipticity of the total mass distribution and the 
modeling, suggesting that the major axis of the potential 
and the external shear differ by PAq — PA c ^t — 20l^g 
deg, that is they are aligned within ~ 1.2a. The orienta- 
tion of external shear in agreement with the orientation 
of stars out to r < 1" which is about —36 deg. The ori- 
entation of the internal quadruple (lens ellipticity) and 
that of stars are misaligned by ~ 24°. Likewise, the axis 
ratio of the light distribution over this radial range is 
0.85 < b/a < 0.93, again consistent with our lens model. 

The lens modeling also puts interesting constraints on 
the redshift of source 2: z S 2 — 5.3 ± 1.0. The accu- 
racy is relatively low because of the saturation of the 
Dig/ D os (z s ) curve when z s -— » oo. The top panel of Fig. 
4 shows a mild correlation between z S 2 and the slope of 
the density profile 7'. This is expected since the steeper 
the density profile that fits the inner ring, the less mass 
is enclosed between the two rings, and hence the further 
away must be the outer source. 

In spite of the complexity of the azimuthal properties 
of the lens potential, our modeling yielded stable and 




Fig. 4.— Top panel: 68.3%, 95.4% and 99.3% CL con- 
tours for model parameters slope of the density profile 7' 
and source 2 redshift z B 2. Bottom panel: Idem for the 
slope 7' and the lens equivalent velocity dispersion (defined 

as lSG^^booq 1 ^ 2 /l" kms" 1 ). 

well localized constraints on the normalization and slope 
of the radial total density profile. The lower panel of 
Fig. 4 shows the confidence regions for the slope 7' and 
the equivalent velocity dispersion asm- First, we find a 
total density profile very close to isothermal with a slope 
7' = 2.00 ±0.03. The corresponding SIE velocity disper- 
sion is cr sie — 287.0 ± 5.2kms~ . In order to compare 
these results with SDSS spectroscopy, one needs to solve 
the spherical Jeans equation taking into account observa- 
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tional effects (SDSS fiber aperture, seeing) and the sur- 
face density of dynamical tracers (radial distribution of 
stars in the lens galaxy) measured in §2.2. Here we as- 
sume an isotropic pressure tensor. A general description 
of the method can be found in Koopmans (2006). Fig. 
5 shows the aperture velocity dispersion that would be 
measured with SDSS spectroscopic fibers when the den- 
sity profile is normalized to fit the first ring alone. It 
shows that slopes close to isothermal (7' ~ 2) predict 
velocity dispersions close to SDSS spectroscopic veloc- 
ity dispersion, which gives strong support to our double 
source plane lensing-only analysis. Such a similarity is 
consistent with the results of previous SLACS studies 
(Treu et al. 2006; Koopmans et al. 2006). We note that 
the accuracy reached on both the slope and the veloc- 
ity dispersion based on lensing constraints alone is bet- 
ter than that afforded by kincmatical measurements at 
the same redshift, although the two methods are comple- 
mentary in their systematic errors and degeneracies (see 
discussion in e.g. Treu & Koopmans 2002). 
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Fig. 5. — Predicted stellar aperture velocity dispersion <r ap 
as it would be measured with SDSS spectroscopic settings as 
a function of the slope of the density profile. The normal- 
ization of density profile is fixed to be consistent with the 
Einstein radius of Ring f . The shaded area is the la SDSS 
measurement uncertainty. It shows a remarkable agreement 
between the double source plane analysis and the coupling 
of kinematical + source I plane data, both favoring nearly 
isothermal slopes. Note that a ap and ctsie do not need to be 
identical. 

4.1. Budget of mass and light in SDSSJ0946+1006 

The tight constraints on the projected mass profile 
between the two Einstein radii can be compared to 
the light distribution inferred in §2.2. In particular, 
the total projected V band mass-to-light ratio within 
the effective radius R c s ~ 7.29 hjQ kpc is M/Ly — 
11.54 ± 0.51 hro (M/Lv) Q (corresponding to a total pro- 
jected mass 4.90 ± 0.13 x 10 11 h^Mg). The logarith- 
mic slope of the projected enclosed total mass profile 
is 3 — 7' = 1.00 ± 0.03, while the slope of the cumu- 
lative luminosity profile close to the effective radius is 
dlogL(< r)/dlogr = 0.62 with much smaller uncer- 
tainty. Therefore the projected mass-to-light ratio profile 
increases with radius as 7-°- 38±0 03 around R c s with high 
statistical significance. 



We now compare these values to the stellar mass con- 
tent in the effective radius using a the typical mass- 
to-light ratio of stellar populations in massive galax- 
ies at that redshift M*/L v ~ 3.14 ± 0.32 h 70 (M/L V ) Q 
(Gavazzi et al. 2007) and ~ 30% intrinsic scatter about 
this values (due to e.g. age-metallicity effects) as found 
in the local Universe (Gerhard et al. 2001; Trujillo et al. 
2004). This leads to a fraction of projected mass in the 
form of dark matter within the effective radius /dm,2d(< 
Reg) — 73 ±9% which is about twice as high as the aver- 
age value found by Gavazzi et al. (2007) and Koopmans 
et al. (2006) thus making SDSSJ0946+1006 a particu- 
larly dark-matter-rich system. 

5. EXPLOITING THE DOUBLE SOURCE PLANE: BEYOND 
THE LENS MASS PROPERTIES 

111 this section we address two particular applica- 
tions afforded by the double source plane nature of 
SDSSJ0946+1006. First, in §5.1 we discuss whether this 
particular system gives interesting constraints on cosmo- 
logical parameters. Then, in §5.2, we present a com- 
pound double lens plane mass model and use it to con- 
strain the total mass of the Ring 1 . This provides a new 
(and perhaps unique) way to obtain total masses of such 
compact and faint objects. Thus, in combination with 
the magnifying power of the main lens, this application 
appears to be a promising way to shed light on the nature 
of faint blue compact galaxies (e.g. Marshall et al. 2007). 
In § 5.3 we discuss the prospects of doing cosmography 
with samples of double source plane lenses, taking into 
account the lensing effects of the inner ring on the outer 
ring. 

5.1. An ideal optical bench for cosmography? 

Can a double source plane lens be used to constrain 
global cosmological parameters like Q m or Q\? In prin- 
ciple this can be done because lensing efficiency depends 
on the ratio of angular diameter distances to the source 
D os and between the lens and the source D\ s as well as the 
projected surface mass density £(0) in the lens plane. In 
formulae, writing the lens potential experienced by light 
rays coming from a source plane as redshift z s as: 

^9, Zs ) = ^-^h f d 2 e'E(o')in\e-e'\ (2) 

c z D OB J 

=M0)^r, (3) 

and considering two images at positions 9\ and 62 
coming from source planes at redshift z s i and z S 2, 
one can measure the ratio of distance ratios r\ = 
(D\ a /D os ) Zs2 /{D\ s /D os ) Zsl directly from the properties of 
the multiple images, given assumptions on the potential 
ipo{&) and its derivatives defining the deflection, conver- 
gence and shear at the positions of the images. 

Applications of this method to clusters of galaxies with 
several multiply imaged systems at different source red- 
shifts - assuming simple parametric models for the clus- 
ters - seem to favor Q m < 0.5 cosmologies (Golse et al. 
2002; Soucail et al. 2004). However, unknown systemat- 
ics lurk under the cluster substructure, which can intro- 
duce significant local perturbations of i/>o(0). In principle 
- at least judging qualitatively from the smoothness of 
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the isophotes, and the smoothness of galaxy scale Ein- 
stein Rings - one could hope that massive elliptical galax- 
ies be less prone to this sort of systematic because source 
size is large compared to the substructure angular scale. 

In the previous section we constrained z S 2 for the given 
ACDM concordance cosmology. Here we re-parametrize 
the problem using 77 itself as a free parameter to con- 
strain the change in lcnsing efficiency between the two 
source planes. The left panel of Fig. 6 shows the joint 
constraints on the two parameters 7' and 77. A first con- 
sequence of this more general parameterization is that, 
by allowing a broader range of values for 77 (i.e. allow- 
ing more freedom in the cosmological model), the un- 
certainties on the slope are significantly increased: we 
find 7' = 2.07 ± 0.06. Steeper density profiles are now 
somewhat compensated by a relatively higher lensing ef- 
ficiency for the second source plane. In other words, the 
tight constraints previously obtained on the slope of the 
density profile depend to some extent on the assumed 
cosmological model (i.e. assuming ACDM cosmology led 
to 7' = 2.00 ±0.03). 
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Fig. 6. — Left panel: Constraints on the logarithmic slope 
7' and the ratio of distance ratios 77. Contours enclose 68.3% 
and 95.4% of probability. The right panel shows r/(z S 2) as a 
function of z S 2 for two flat cosmologies (fi m ,fiA) = (0.3,0.7) 
(black) and (H m ,fiA) = (l.,0.) (red) which are two sensible 
"extreme" cases. The dotted horizontal lines illustrate the 
upper limits on 77 for these cosmologies given the assumption 
z s2 < 6.9 (see §2.3). 

The right panel of Fig. 6 shows 77(2,52) as a function 
of the second source redshift for two "extreme" flat cos- 
mologies: (n m ,n A ) = (0.3,0.7) and (O m ,0 A ) = (l.,0.), 
intermediate cases lying in between. This shows that 
high values 77 > 1.57 are not consistent with currently 
favored cosmologies. The upper limit on ?7(z S 2 = 6.9) is 
also shown for these two cases. This illustrates that very 
loose constraints can be obtained on cosmological param- 
eters even if z S 2 were known spectroscopically. Likewise, 
even assuming an isothermal slope of the density pro- 
file as motivated by joint lensing and dynamical analyses 
(Koopmans et al. 2006) does not drastically improve the 
constraints on 77 and consequently on cosmology as shown 
in Fig 7, even if z s2 could be measured with spectroscopic 
precision. 

However, it is important to point out that the formal 
~ 3% relative uncertainty we get on 77 from our lens 
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Fig. 7.— 68.3%, 95.4% and 99.3% CL contours in the red- 
shift of source 2 and Q m parameter space assuming an isother- 
mal density profile. This shows that even using strong priors 
on the density profile and for a given source redshift, only 
loose constraints can be inferred on cosmological parameters 
with a single double source plane system. 

modeling strategy based on the identification of conju- 
gate knots underestimates the potential accuracy of the 
method. Statistical errors would decrease by a factor 
of a few with a full modeling of the surface brightness 
distribution in the image plane. Unfortunately, the er- 
ror budget would then be limited by additional system- 
atic sources of uncertainty like extra convergence coming 
from large scale structures along the line of sight with 
estimated standard deviation er K > 0.02 (Dalai et al. 
2005) or due to a non trivial environment in the main 
lens plane. Therefore, we conclude that it is unlikely 
that any cosmographic test based on the unique multi- 
ple Einstein ring system SDSSJ0946+1006 will provide 
valuable information on cosmological parameters. The 
prospects of using large numbers of double source plane 
lenses are investigated in §5.3. 

5.2. Source 1, alias Lens 2 

Among the massive perturbers along the line of sight 
to source 2, the most prominent is probably the mass as- 
sociated with source 1. Since the lens modeling predicts 
that both sources are located very close to the optical 
axis (the center of the lens, see lower panels of Fig. 3), 
the light rays coming from source 2 to the observer will 
experience the potential of source 1 before that of the 
main lens. Fig. 8 illustrates the complexity of the con- 
figuration which adds some extra positive focusing for 
the second source plane. For the cosmological applica- 
tions we mentioned above, this translates into a small 
but systematic source of bias. The bias introduced on 
the inferred mass profile of the main lens is small, so 
that the conclusions presented in §4 are not significantly 
altered except for the estimate of z S 2- 

On the bright side, this lens configuration allows us to 
obtain some insight on the mass associated with Ring 
1 (also identified as "Lens 2") provided we now fully 
take into account the multiple lens plane nature of such 
lines of sight (e.g. Blandford & Narayan 1986; Schneider 
et al. 1992; Bartelmann 2003). This is the purpose of the 
present section, in which we fix the ACDM concordance 
cosmological model for simplicity. 
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Fig. 8. — Sketch of the lensing optical bench with source 
1 acting as a perturbing lens on source 2 which complicates 
the relation between redshifts, deflection angles and angular 
distances. 

To achieve this goal we have to address the mass prop- 
erties of the main lens at the same time as those of the 
first source 1. We reconsider the lens model of §3, but add 
another mass component at redshift zn = z s i = 0.609 in 
the form of a singular isothermal sphere with free equiv- 
alent velocity dispersion parameter asiE.si and centered 
on the position of source 1. As in §3, our uncertainty 
on the distance to source 2 is simply parameterized by 
its redshift z S 2 in the context of a ACDM cosmological 
model. 

In multiple lens-plane theory, the relation between the 
angular position 6j of a light ray in the j-th lens plane 
and the angular position in the j = 1 image plane is: 



i = l 



■•)• 



(4) 



The last lens plane N can be identified with the source 
plane such that 9m = f3. In Eq. 4, as compared to 
Bartelmann (2003), we did not consider the reduced de- 
flection which introduces an unnecessary extra Di S /D s 
term in the sum. Likewise, the sign convention for the 
deflection is different than Bartelmann (2003). There- 
fore for two distinct positions 9\ and 02 coming from two 
distinct source plane positions f3± and j3% respectively, we 
can write: 
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(5) 



In these equations, a is the deflection produced by the 
main lensing galaxy (lying in the plane that also defines 
the image plane) and a s i is the perturbing deflection 
produced by source 1 (lens 2) onto source 2. Note that 
parameters like the center of source 1 enter the modeling 
scheme both as source- and lens-plane parameters. This 
is clearly visible in the brackets for the argument of a s i 
that contains /3i, the position of source 1 in the source 
plane. 

The constraints obtained on the equivalent velocity 
dispersion parameter of the main lens ctsie and that 
of source 1 ctsie, s i are shown in the left panel of Fig. 
9. We clearly see two kinds of solutions: one (fam- 
ily i) has a high lens velocity dispersion (and slope 
7' ~ 1.96, nearly isothermal) and little mass in source 
1, whereas the other family (ii) has a lower main lens ve- 
locity dispersion and more mass in source 1. We measure 
(VsiE) csiE.si) = (295^5 0, 56±30) kms -1 for family i and 
(oSiE.ffsiE.ai) = (247.3t^,104;2i)kms- 1 for family ii 



2.3 2.2 2.1 2.0 1.9 1. 

I " " I " " I " " I " " I " 



150 - 



100 



50 




T"T"T"T"T 




150 



100 



50 



220 240 260 280 300 32Q 



-'SIS.lens 



[km/s] 



source 2 redshift z s2 



Fig. 9. — Left panel: contours in parameter space of 
the velocity dispersion of the main lens ctsie and that of 
the first source ctsie, si. Given the tight correlation ctsie — 
(687 — 200.37') kms -1 found in §4, the upper abscissa shows 
the correspondance with slope 7'. The kinematical SDSS es- 
timate of CTu,» and the velocity dispersion of source 1 inferred 
from the Tully-Fisher relation (Moran et al. 2007) are over- 
laid as a point with error bar. Right panel: contours in 
parameter space of the second source redshift z B 2 and the ve- 
locity dispersion of the first source ctsie, si- The recovered 
z S 2 strongly depends on the mass enclosed in source 1. In 
both panels confidence levels mark the 68.3, 95.4 and 99.3% 
enclosed probability. 

A pixelised source plane inversion for both of these best 
fit models is shown in Fig. 10. Family i models are shown 
in the top row and family ii in the bottom row. Note the 
very complex systems of caustic and critical lines pro- 
duced by this multiple lens plane system. It is difficult 
to favor either of these models based on a visual inspec- 
tion and either region on the parameter space has about 
the same statistical weight (fraction of MCMC samples). 




Fig. 10. — Top panels: Best fit family i model image and 
source plane reconstructions. From left to right reconstructed 
image plane, residual (data-model), and source 2 plane at 
redshift z S 2 = 3.30. Bottom panels: idem for the best fit fam- 
ily ii models (with z S 2 = 2.75). Note the complex critical 
and caustic curves for the z S 2 source plane due to the multi- 
ple lens plane configuration produced by source 1. For both 
models the reconstruction is satisfying and produces very few 
residuals. 
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The left panel of Fig. 9 also shows the aperture- 
corrected SDSS-inferred velocity dispersion of the lens 
cr„ ; * = 284 ± 24kms~ 1 which seems to favor family i 
solutions, based on the earlier SLACS results of a gen- 
eral agreement between stellar velocity dispersion and 
asm- In addition, we can get further external infor- 
mation on the mass of source 1, by extrapolating the 
Tully-Fisher relation found by Moran et al. (2007) at 
z <~ 0.5 for late- type galaxies. In the field, they found 
that at absolute magnitudes of V ~ —19.7, the maxi- 
mum rotation velocity is log(2y max ) = 2.2 ± 0.1. As- 
suming V^nax ~ v^ctsie, this translates into an estimate 
<7SiE,si — 59 ± 13 kms -1 . Another piece of information 
comes from weak lensing results at intermediate redshift 
(0.2 < z < 0.4) by Hoekstra et al. (2005), who found 
that galaxies with magnitude V — 51og/i ~ — 19 have 
virial masses M vir ~ 1.50+°;g4 x 10 11 ftf 1 M Q which also 
corresponds to log(2V max ) = 2.20 ± 0.09, in good agree- 
ment with Moran et al. (2007). These two arguments 
also seem to favor family i solutions, i.e. those with more 
mass in the main lens and less in source 1. 

The right panel of Fig. 9 shows the important degen- 
eracy between the redshift of source 2 and the velocity 
dispersion of source 1. We can see that the more massive 
source 1, the lower z S 2 must be. This demonstrates that 
any cosmographic test based on multiple source plane 
lens systems should carefully consider the mass in the 
foreground source as a significant perturbation on light 
rays coming from the most distant source. Adding a sub- 
stantial amount of mass in source 1 significantly changes 
the inferred redshift of source 2 either for family i models 
which yields z s2 = 2.6+J;° or family ii models yielding 

z s2 = 3.81^5 . Marginalizing over the whole posterior 

PDF gives z s2 = 3.1±?;g . 



TABLE 3 

Best-fit model parameters for SDSSJ0946+1006 using a 

COMPOUND DOUBLE LENS PLANE. 
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Best fit model parameters and 68.4% confidence limits. An- 
gles are in degrees oriented from North to East. 

5.3. Future outlook: cosmography with many double 
source plane lenses 

In §5.1 we explored the possibility of constraining cos- 
mology with SDSSJ0946+1006, and came to the conclu- 
sion that the errors are too large for this to be interest- 
ing. In §5.2 we saw that the mass of the closest source 
must be taken into account as a perturbation along the 
double source plane optical bench. Here we attempt to 
address the possibility of using large numbers of such 



multiple lensing systems to probe the cosmology. Future 
space-based missions like DUNE or JDEM should pro- 
vide us with tens of thousands of lenses, among which 
several tens would be double source plane systems. We 
also assume that redshifts will be available, from space- 
or ground-based spectroscopic follow-up. 

First, we summarize the error budget expected for a 
typical double source plane system. As described before, 
the main quantity of interest is the ratio of distance ratios 
parameter 77 = {D\ s / D os ) 2 / {D\ s / D os )i, where source 2 is 
the furthest one. For simplicity, we assume that the main 
lens, the first source, and the second source are perfectly 
aligned onto the optical axis, resulting in two complete 
concentric rings of radius Q\ and 6 2 . The lens equation 
for each source plane reads: 

(3i=e 1 -(D ls /D os ) 1 a tot {6 1 ) = 0, (7a) 
02 = 62 -{Du/D„) 2 a tot (02) = Q. (7b) 

We consider again the general power-law surface mass 
distribution of Eq. (1) producing deflections ct\ and a 2 
on source 1 and source 2 light rays. For source 2 we must 
add a p the small perturbing deflection 11 due to source 1 
and experienced by source 2 only. Combining Eq. (7a) 
and Eq. (7b) gives: 

" = (I) i + a^an • (8) 

This equation shows the importance of the perturba- 
tion. If one aims at constraining 77 with interesting accu- 
racy (i.e. error smaller than 0.01), the small perturbing 
term in the denominator of second part on the right hand 
side of Eq. (8) should be smaller than 0.01. Keeping 
in mind that for lensing potentials close to isothermal, 
a oc a 2 , and that the typical velocity dispersion of the 
main lens is about a ~ 250kms _1 , it is important to 
control and correct perturbing potentials with velocity 
dispersion as small as cr p — aj 10 <~ SOkms^ 1 for values 
D sls2 /D ls2 ~ 0.5. 

Next, differentiating Eq. (8), and writing r = 9 2 /9i, 
one can infer the fractional error on rj: 

(^) 2 = (7'-D 2 (^) 2 + (ln^ + 



1/ + D sla2 ^ p ) 

The first contribution is the relative measurement er- 
ror on the ratio of Einstein radii, with typical values 
0.001 < S r /r < 0.03 for deep space based imaging. The 
second term captures our prior uncertainty on the slope 
of the density profile (for example Koopmans et al. (2006) 
measured (7') ~ 2.01 and an intrinsic scatter 6 7 > ~ 0.12). 
Finally, the third term represents our prior knowledge of 
the mass of the perturber, which can be based, for ex- 
ample, on the Tully-Fisher relation. Moran et al. (2007) 

11 We assume that the non-linear coupling between lens planes 
can be neglected, i.e. the perturbation of source 1 is small com- 
pared to the deflection from the main lens on source 2 light rays: 
ctp <C 012 ~ O2 ■ 
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estimated S ap /<r p ~ 0.25. Inserting those values into (9), 
and assuming a typical value of r ~ 1.5, we find a rela- 
tive uncertainty on S^/rj ~ 0.06 for a single system. The 
error is dominated by model uncertainties on the slope 
of the density profile and the mass in source 1. In the 
case of SDSSJ0946+1006, we achieve a similar accuracy 
when using the above priors on the slope and the veloc- 
ity dispersion of source 1. In the following we shall use 
a conservative S^/r) = 0.08 fiducial value. 

Having estimated the accuracy achievable on r\ for a 
single double source plane system, we focus on the cosmo- 
logical meaning of rj in a spatially flat universe dominated 
by Dark Matter (O m ), and Dark Energy (Ode = 1 — ^m) 
with equation of state parameter w = pde/pde- It is 
worth pointing out that the ratio of angular diameter 
distances is independent of the Hubble constant Hq. 

The constraints on (Q m ,w) are shown in the right 
panel of Fig. 11. The error contours are obtained us- 
ing a Fisher matrix formalism. We assumed a sample 
of 50 double source plane lenses, randomly produced us- 
ing Monte-Carlo simulations. The redshift distribution 
of lenses and sources used for the Monte-Carlo simula- 
tions is shown in the upper left panel of Fig. 11. For 
the parent population of sources, it is based on recent 
COSMOS estimates (Leauthaud et al. 2007). The equiv- 
alent velocity dispersion of the lenses is assumed to be 
Gaussian with mean and standard deviation of 190 and 
60kms _1 respectively. The Einstein radii for the first 
and second ring are constrained be greater than 0.7" and 
1.0" respectively. 




Fig. 11. — Top left panel: Redshift distribution for the lens, 
source 1 and source 2 used in the Monte-Carlo simulation of 
mock double source plane lenses as they could be discovered 
in future space-based surveys. Bottom left panel: Individual 
error brought by each system on Sl m as a function of the 
redshift of the first source. Right panel: 68.3, 95.4 and 
99.3% confidence levels contours on the matter density f2 m 
and equation of state of dark energy w = p/p obtained when 
combining 50 multiple source plane systems. 

The cosmological parameters Cl m and w are recovered 
with a precision ±0.020 and ±0.080 respectively. We 
note that the sensivity and the orientation of the degen- 
eracy in this set of cosmological parameters is similar 
to those obtained with a type la supernovae experiment 
(see e.g. Refregier et al. 2006, and references therein). 
The lower left panel of this figure demonstrates that the 



systems that contribute the most to constraining f2 m are 
those with source 1 redshift z s i > 1 (a similar trend is 
seen for w). Lens redshifts larger than ~ 0.5 are also 
more efficient configurations. This can easily be under- 
stood since the higher the redshifts of either lens, first or 
second source the more sensitive distances are on cosmo- 
logical parameters. Note that situations with very low 
redshift lenses but high redshift source 1 and 2 will re- 
sult in a rapid saturation of the D\ s /D os for both source 
1 and 2, leading to values 7/ ~ 1 independant on cosmol- 
ogy. The sensitivity on cosmology is actually essentially 
driven by the redshift of the primary lens. Therefore we 
conclude that SDSSJ0946+1006 is not an optimal double 
source plane lens system for cosmographic purposes. 

However the prospects of doing cosmography with fu- 
ture samples of double source plane lenses are excel- 
lent, provided systematic effects are controlled. The 
main source of systematic uncertainty that was ignored 
in the above calculations is the possibility of a change 
in the mean density slope as a function of lens red- 
shift (7') = f(zi) as present data seem to suggest that 
the dynamical properties of early-type galaxies have not 
changed much since z ~ 1. Koopmans et al. (2006) found 
d( 7 ')/dz = 0.23±0.16over the redshift range [0.08,1.01]. 
Progress needs to be made along this line to improve how 
knowledge of (7') = f{z\), but a great advantage of dou- 
ble source plane systems over single ones is that com- 
bining stellar dynamics and lcnsing constraints from two 
source planes would be more efficient at "self-calibrating" 
the method than using single ones. In addition a thor- 
ough lensing analysis aiming at carefully modeling the 
surface brightness of lensed structures will certainly help 
in controlling any such evolution trend of the density 
profile (see e.g. Dye & Warren 2007). 

6. SUMMARY & CONCLUSION 

In this paper we report the discovery of the first galaxy- 
scale double lensing event made of a foreground lens 
galaxy at redshift z\ = 0.222, a first source at redshift 
z s i = 0.609 (Ring 1) and a more distant source (Ring 2) 
with unknown redshift, despite an attempt to measure 
its redshift with deep optical spectroscopy using LRIS 
on the Keck I Telescope. The detection of Ring 2 in a 
single orbit HST-ACS F814W filter image, sets an upper 
limit to its redshift z S 2 < 6.9. 

Modeling the geometry of the lensed features at differ- 
ent source planes we determine the mass density profile 
of the lens galaxy which is found to be close to isother- 
mal. The best fit lens model predicts a stellar veloc- 
ity dispersion in very good agreement with that mea- 
sured from SDSS spectroscopy. The model requires a 
relatively large amount of dark matter inside the ef- 
fective radius /dm,2d(< Res) — 73 ± 9% (correspond- 
ing to a projected total mass-to-light ratio M/Ly — 
11.54±0.13 hro (M/Ly) Q ), assuming the stellar mass to 
light ratio measured in paper IV. Along with the com- 
plex isophotes of the lens galaxy and the presence of 
several other (less luminous) galaxies at similar photo- 
metric redshifts, the high dark matter fraction suggests 
that the lens may be the central galaxy of a group scale 
halo. The high precision of this measurement - far su- 
perior to that attainable from a single multiply imaged 
systems - demonstrates that double source plane lenses 
are extremely valuable tools to study the mass profile of 
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galaxies and groups. 

In order to constrain the redshift of Ring 2 and as- 
sess the feasibility of determining cosmological param- 
eters using double source plane lenses, we constructed 
multiple lens plane mass models. In that case the lens- 
ing effect of Ring 1 on Ring 2 is taken into account and 
modeled as a singular isothermal sphere. Although the 
extra mass component adds additional uncertainty to the 
derived z s2 and cosmological parameters, it provides a 
unique way to determine the total mass of the interme- 
diate galaxy (the inner Ring). Discarding the family ii 
range of solutions (disfavored by kinematics of stars on 
the main lens galaxy), the two- lens-plane mass model 
results can be summarized as follows: 

• The redshift of Ring 2 is found to be z s2 — 2.6+ ' 7 . 
This is a genuine prediction that can be tested with 
the help of deep Hubble images at shorter wave- 
lengths, using the dropout technique. 

• No interesting constraints on cosmological param- 
eters can be obtained from the lensing analysis of 
the system SDSSJ0946+1006, due to the unknown 
redshift of Ring 2, the overall degeneracy of cos- 
mography with the slope of the mass density profile 
between the rings, the degeneracy with the mass of 
the inner ring, and the suboptimal combination of 
lens and source redshifts. 

• The velocity dispersion of Ring 1 is found to be 
csie = 56±30kms~ 1 , in good agreement with the 
value expected based on the extrapolation of Tully- 
Fishcr relation at this redshift (Moran et al. 2007) 
and on weak-lensing measurements (Hoekstra et al. 
2005). Given that lensing cross section increases 
with the fourth power of the velocity dispersion, 
individual lenses in this mass range are expected 
to be very rare. In addition, it is observation- 
ally more difficult to identify lensed background 
structures embedded in foreground late- type galax- 
ies with small Einstein radius (both for imaging 
or spectroscopic lens searches). Thus by exploit- 
ing the boost of the primary lens, double source 
plane, also seen as double lens plane systems, may 
be an effective way to determine the lensing mass 
of small distant galaxies, complementing detailed 
photometric studies (Marshall et al. 2007), and 
kinematic studies with integral field spectrographs 
on large ground based telescopes with adaptive op- 
tics. 

Future planned space missions like JDEM or DUNE 
are expected to deliver several tens of thousands of single 
source plane lenses (Aldering & the SNAP collaboration 
2004; Marshall et al. 2005; Refregier et al. 2006) and 
several tens of double source plane lens galaxies. Given 
the great utility of multiple source plane lenses as tools to 
study distant galaxies, and the relatively small number 
of expected systems, we argue that the necessary effort 
of spectroscopic follow-up would be easily affordable and 
well motivated. 

In addition, a relatively large sample of double source 
plane galaxy-scale gravitational lenses will be a practical 



tool for cosmography. As an example, we calculated the 
constraints on £l m and the equation of state of Dark En- 
ergy w = pde / Pde that can be obtained from a sample of 
50 double source plane lenses, assuming both source red- 
shifts are known and are realistically distributed. Spec- 
troscopic follow-up of such systems rings is also required 
to control systematic effects such as the change in the 
mean density profile slope as a function of the lens galaxy 
redshift. A careful analysis taking into account the un- 
certainty on the mass profile of the main lens and of 
the perturber shows that cosmological parameters can 
be measured with an accuracy of 10% comparable to 
that obtained from the Hubble diagram of Type la su- 
pernovae. 
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APPENDIX 

DERIVATION OF THE PROBABILITY OF MULTIPLE LENSING 

In this Appendix we estimate the probability of finding a double lens in a sample of lenses like SLACS. The first 
ingredient is the surface density on the sky of potential lens galaxies, given by: 

A^/d./d^O^, (Al) 

with ^ the comoving volume per unit solid angle and redshift, and the velocity dispersion function. For 
simplicity, we assume here that the shape of the velocity dispersion function does not evolve with redshift, but only in 
normalization as described by the p(z\) function. In practice, we consider the velocity dispersion function ^ measured 
by Sheth et al. (2003) at z <~ 0.1, which is of the form: 

£^)°*m->W'tf- (A2) 

with 0* = 0.0020 ±0.0001 hf Mpc~ 3 , cr* = 88.8 ± lT.Tkms" 1 , a = 6.5 ± 1.0 and (3 = 1.93 ±0.22 . 

The number density of foreground galaxies producing a single strong lensing event on a source population si can be 
written as: 

f f f dV dn; / \diV s i 
N sl = J dz s i J da J dzi— —p(zi)-j—X{a,z h z s i) , (A3) 

following Marshall et al. (2005). In this equation, X(a, z\, z s \) is the cross-section for lensing. In most cases of scale- 
free gravitational lenses X(a,zi,z s i) can be separated such that X(a, z\,z s \) = a 2 " g{z\, z s \), a v giving the overall 
strength of the lens. For the particular case of a singular isothermal sphere that we shall consider, v — 2 and 
g oc (Di s i/£> os i) 2 0(z s i — z\), with Q(x) the Heaviside step function. 

If the lensing cross-section for a second population of sources s2 does not depend on the presence or properties of 
an already lensed population si galaxy, we can write: 

Nsl, s2 = j dz s2 J dz s i J da J dzi^^p{z x )'^^'^^X l {a, z h z sX )X 2 {a, z x , z s2 ) . (A4) 

Combining Eqq. (Al), (A3), (A4) and taking advantage of the separability of the dependency on a and on redshifts, 
the ratio of the probability that a galaxy lenses a source at z s2 given that it is already lensing a source at z s \ over the 
probability for a galaxy to lens a source at z s2 is given by: 



(A5) 
(A6) 



P(lenS S2 | lens Si) _ P(lens s2, lens si) _ Nsl^iVga! 

P(lcnS S2) ~~ P(lens s2) P(lons si) — N a2 N s2 

_ [f da^] X \f da^a^] [f dVjzQp^)] X [j d Zal ^L J d Zs2 ^ J dVjztMzQgiz^Mz^)] 

~ [fda^-^Y [/ dz 3l ^ / dV^MzMzuZ.l)] X [/ dz B2 ^f f dV( Zi )p(zi)s(zi,z B 2)J 

= Ex(. (A7) 

where dV(zi) indicates ^dzi. The first term S in Eq. (A6) describes the strong dependency of the lensing cross 
section on velocity dispersion. As expected because lensing favors high a systems, using the velocity dispersion function 
from Sheth et al. (2003), we estimate it to be larger than unity, of order S = r((8 + a)/ / 3)r(a//3)/r((4 + a)//3) 2 ~ 2.44. 

The second term £ contains volume and lensing efficiency g{z\, z s ) effects that depend on the redshifts of the lens 
and the sources. By defining a lensing efficiency averaged over a given population of sources, 

Gi(zi) = j &z si ^-g{z\,z si ) , (A8) 

we can simplify the second term in Eq. (A6) and write it as: 

c = [JW(f0pW] x [/dV(z ; )p(z 1 )G 1 (z 1 )G 2 (z 1 )] 

[/dV(*,)p(*i)Gi(zi)] x [dV{zi)p{z^G 2 {z{)\ ' ( ' 

To obtain a quantitative estimate of the probability of double lensing, let us consider some specific examples. The 
most important quantity in the definition of £ is the comoving redshift distribution of deflectors p(z\). If all of them 
were confined in a single lens plane such that p(z\) = 5{z\ — z\q), then Q = 1 and most of the change in probability 
comes from selection effects captured by E. If, instead, deflectors are broadly distributed over a range of redshifts, 
the ratio can be significantly higher than one because the probability of single lensing would be low for high redshift 
deflectors whereas the fact that a deflector is already lensing a source at z s \ favors lenses in a low redshift range, more 
suitable for lensing a source at z s2 . 

We illustrate this volume effect by assuming that the comoving density of deflectors is constant out to a redshift z max 
and then it drops to zero, that is p{z\) = 9(2; ma x — z\). We also assume the redshift distribution of background sources 
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is of the form cUV/dz oc e~ z / z °(z/zo) a ~ 1 . Population s2 galaxies follow the redshift distribution of faint background 
sources presented in Gavazzi et al. (2007) and having zq = 0.345 and a — 3.89. We use a redshift distribution for the 
si population that peaks around redshift 0.5, in agreement with the properties of spectroscopically discovered SLACS 
lenses (see Bolton et al. 2008). This corresponds to z ~ 0.07 and a ~ 7. Note that the detailed shape of the redshift 
distribution for either si or s2 galaxies does not change the trends significantly, so this approximation is sufficient 
for our purposes. Fig. Al shows the evolution of ( as a function of the limiting redshift z max . We see that, out to 
reasonable values z max ~ 0.5, ( does not depart much from unity. 

To obtain a numerical value to be compared with our SLACS sample, we consider a population of deflectors constant 
out to redshift unity. The gain in probability is P(lens s2 | lens sl)/P(lens s2) ~ 2.4 — 5. In other words, if one elliptical 
galaxy at z < 0.8 in about 200 is strongly lensing a faint background source, a strong lens in approximately 40-80 is 
a double lens. This is consistent with the observations. 




Fig. Al. — Evolution of the £ term for multiple lensing probability boost P(lens s2 | lens sl)/P(lens s2) as a function of the 
limiting redshift 2 max of the lens distribution. 
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